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Abstract 

An  infinite-dimensional  nodal  la  given  foe 
the  ganaratlon  of  gyroscopic  noisa,  which  exhibits 
power  spectral  density  proportional  to  (1/f)  over 
a wide  frequency  range.  The  optimal  filter  is 
given  for  separating  a statistically  described  sig- 
nal froa  additive  gyroscopic  noisa,  using  dlscrate- 
tiae  observations.  This  filter  is  oppressed  as  a 
discrete- tine  infinite-diaensional  KAlnan-Bucy 
filter,  with  an  associated  Riccati  covariance 
operator  equation.  Sufficient  conditions  era  spe- 
cified such  that  this  Kalaan-Bucy  filter  will 
possess  various  desired  properties. 


1.  Introduction 

The  gyroscope  is  an  lnstnaant  uaed  to  detect 
angular  movement.  The  problaa  of  the  reaoval  froa 
the  gyro  output  signal  of  noise  Inherent  to  the 
gyroscope  in  a constant  gravitational  field  is  one 
which  has  received  considerable  attention  in  the 
literature.  Sutherland  and  Gelb  tl] , for  exaaple, 
discuss  an  aided  inertial  guidance  system,  where 
periodic  telescopic  sightings  are  uaed  along  with 
gyro  output  to  develop  gyro  error  observations. 

The  error  observations  are  used  as  the  input  to  a 
Kalman  filter,  which  is  used  to  estimate  the  gyro 
error  at  the  observation  times.  An  estimate  of 
the  true  angular  position  is  then  obtained  by  sub- 
tracting tha  estimated  gyro  orror  from  the  gyro 
output  samples.  Mehra  and  Bryson  (2)  discuss 
smoothing  of  the  gyro  output  to  obtain  estimates 
of  the  input  signal. 

Gyroscopic  noisa  has  often  been  modeled  as 
either  a first-order  Causa-Markov  process  (3),  or 
as  a Gaussian  random  walk  (Integral  of  Gaussian 
white  noise)  14,51.  However,  recant  studies  per- 
formed at  The  Charles  Stark  Draper  Laboratory  (6)  * 

of  the  power  spectral  characteristics  of  the  random 
noise  associated  with  various  gyroscopes  indicate 
.hat  gyro  noise  is  often  characterised  by  a (1/f) 
behavior  in  power  spectral  density  over  a wide 
frequency  range.  An  explanation  of  the  source  of 
this  noise  in  the  magnetic  materials  of  the  gyro- 
scope (e.g.  the  gyro  float  rebalance  torquar)  is 


The  research  of  tha  first  author  was  supported  by 
the  National  Science  foundation  under  a Graduate 
Fellowship.  The  research  of  the  second  author 
was  supported  by  the  Air  Force  Office  of  Scientific 
Research  under  grant  AF-AF05R-72-2273. 


proposed  by  Harris  and  Koenigaberg  (7) . h Sec- 
tion 2 we  discuss  their  findings  and  add  others. 

Me  present  an  infinite-dimensional  state  space 
model  which  generates  noise  with  the  power  spec- 
tral properties  of  gyroscopic  noise.  Me  also 
discuss  the  possible  relationship  between  magnetic 
dlsaccn— mdatlon  and  gyroscopic  noise. 

In  Section  3 we  introduce  and  solve  the  fil- 
tering problem  to  be  treated  in  the  paper.  Using 
discrete-time  observations,  a statistically 
described  gyro  output  signal  (resulting  from  angu- 
lar motion  inputs  to  the  gyroscope)  is  optimally 
separated  from  additive  gyroscopic  noise.  Be- 
cause observations  are  made  at  discrete  times,  we 
first  determine  a discrete-time  infinite-dimen- 
sional linear  system  to  generate  samples  of  the 
gyroscopic  noiss,  as  modeled  in  continuous  time  in 
Section  2.  The  filtering  problem  can  be  solved  as  a 
conditional  expectation  filter  in  the  case  where 
the  input  signal  is  Gaussian  (this  solution  being 
equivalent  to  the  minimum  variance  linear  estima- 
tor for  non -Gaussian  input  signals) . The  result- 
ing optlmsl  filter  is  expressed  as  a discrete-time 
infinite-dimensional  Kalman  filter  with  an  asso- 
ciated Riccati  covariance  operator  equation.  Me 
note  here  that  steady-state  filtering  of  a random 
process  with  a (f-l-2e>  power  spectrum  has  been 
diseuasad  by  Moran  (8) . However,  the  performance 
of  Moran's  fllear  degrades  as  C ♦ 0. 

Ns  indicate  how  theorems  concerning  Hilbert 
space  Kalman  filters  and  Riccati  operator  equations 
can  be  appliod  to  the  gyro  noise  filtering  problem. 
By  specifying  conditions  on  the  system  generating 
tha  signal  to  be  recovered,  we  are  able  to  guaran- 
tee a mater  of  desirable  properties  for  the  Kal- 
mar. filter. 

The  optimal  filter  derived  in  Section  3 in- 
volves integrations  over  a free  time  constant  para- 
meter. at  applications,  these  integrations  must 
be  irplemented  discretely.  This  discretisation 
can  be  achieved  by  making  a finite-dimensional 
approximation  to  the  infinite-dimensional  gyro- 
scopic noise  model.  The  optimal  filter  bscoetes  an 
ordinary  finite-dimensional  discrete-time  Kalman 
filter,  with  an  aaaociated  matrix  Riccati  equation. 
It  can  be  shown  (19)  that  tha  mean-squared  estimation 
error  Incurred  in  using  the  Kalman  filter  of  the 
finite-dimensional  approximate  model  can  be  made, 
through  the  use  of  a sufficient  number  of  dimen- 
sions in  the  approximation,  to  approach  the  mean- 
squared  estimation  error  associated  with  optimal 
filtering  of  gyrosopic  noise. 


Proceedings  1976  Decision  and  Control  Conference,  CleaiViter  Beach,  Florida, 
December  1976, 


fPF 


‘i'll  M»s  ' t | .•  J • »••*.  ’ . ,1  )• 


f 


Dm  (Malta  found  here  for  gyroscopic  aoiaa 
mo  applicable  to  any  random  process  charactorlcod 
by  a wideband  (l/f)  powar  apactral  danalty.  aa 
Viuon  in  Suction  2.  (l/f)  spectra  ara  found,  for 

•uaapia  in  aaatlconductor  flickar  no  I aa  and  In  tha 
aoiaa  characterising  tha  fra<|uancy  fluctuation  of 
quarts  cryatal  oeclllatore. 

2.  an  fcfinlta-Plnanaional  modal 
for  Oyroacoplc  Moiaa 

Cyroacopic  noiaa  haa  of tan  boon  aodalad  aa 
aithar  a firat-ordor  Gauaa-Harkov  procaaa  (11,  or 
aa  a Gauaaian  random  walk  (Integral  of  Cauaalan 
whita  noiaa)  (4,51.  Tha  Gauaaian  natura  of  tha 
Miaa  ia  infatrad  from  hlatogra*  plota  of  gyro 
output.  A linoarisad  version  on  log-log  acalaa  of 
tha  powar  apactral  danaity  of  a flrot-ordar  Gauas- 
Nnrhov  procaaa  la  ahown  in  Pigura  1.  Tha  random 
walk  hu  a variance  proportional  to  tima,  banco  la 
nooatatlonary.  Thua  in  a atrlct  aanaa  tha  powar 
apactral  danaity  of  a random  walk  procaaa  dooa  not 
aniat.  Xhen  diacrata  aaaplaa  of  band limited  white 
noiaa  ara  generated  by  conputar  and  s«ad  (to 
reenable  tha  integration  of  whita  noiaa) , tha  re- 
acting noiaa  la  found  to  be  characterlsad  by  a 
(l/fl)  powar  apactral  danaity  over  tha  benbaldth 
of  tha  original  bandlimltad  whita  noiaa.  (Tha 
power  apactral  danaity  la  found  through  evaluation 
of  tha  squared  nagnitudea  of  tha  Fourier  coeffi- 
cienta  of  tha  output  aignal.)  For  tha  following 
raaaona  wa  Intuitively  aapact  thla  reault.  Tha 
powar  apactral  danaity.  Syy(f).  of  tha  output  of  a 
time-invariant  linear  ayatam  (tranafer  funcelona 
0(f) ) to  an  input  aignal  of  FSD  (powar  apactral  dan- 
aity) B„(f)  ia  given  byi 


•py <f)  - •**(« 


|«(f)|2 


(1) 


tranafer  f wet  ion  of  aa  Integrator  la  propor- 
to  (1/a),  hanca  wa  would  havei 

. Jk 1 

()2wf)  l-)2»t)  TrV 


(2) 


> will 


Email  ini  tad  whita.  noiaa  haa  a FSD  oonatant  with 
frequency  (over  ita  band  llnita) , ao  we  would  in- 
tuitively expect  our  approximation  to  random  walk 
to  have  behavior  proportional  to  (l/f2) . Tha  FSD 
remelting  from  the  coaputer  simulation  deecrlbed 
Wove  la  vkff-m  In  Figure  2.  notice  that  both  ran- 
dom proceaaM  dlacuaaad  hara  exhibit  (l/f2)  behav- 
ior in  PSD  (slopM  of  (-2)  on  log-log  ocalee) . 

Decent  etediee  performed  at  Tha  Charlee  Stark 
Draper  Laboratory  It)  of  tha  powar  apactral  char ac- 
tor latica  of  tho  random  noiaoa  aaaoelatad  with 
varioua  gytoaoopaa  indicate  that  gyro  noiaa  la 
often  charactarlaad  by  a (l/f)  behavior  in  powar 
apactral  daaalty.  (Tha  gyro  ia  eat  up  m aa  input 
rata  integrator,  with  a binary  torque  rebalance 
loop.  Tho  unite  of  tSD  are  (input  rata) 2/Xa. ) A 
linearised  graph  of  tho  oboorved  form  of  gyro  power 
spectral  density  is  given  in  Figure  ).  The  (f2) 
portion  of  this  tfraph  ia  primarily  attributed  to 
grant  laatioo  noiaa  due  to  tho  binary  torque  1«|. 
This  affect  of  quantisation  is  currently  under  in- 
vestigation. tower  apactral  analyaaa  of  asperate 
record  langthe  of  gyro  noiaa  show  tho  power  ap 
tral  danaity  to  bo  oonatant  to  time,  hence  wa 


of  tho  eource  of  thla  noiaa  in  tha  magnetic  mater- 
lala  of  tha  gyroscope  (e.g.  tha  gyro  float  reba- 
lance torquer)  ia  propoaad  by  Harris  and  Koenigs- 
berg  (71 . In  thla  section  wa  shall  discuss  their 
findings  and  add  others.  Ha  first  discuss  a modal 
for  magnetic  relaxation- (dlaaccommodatlon) . Thla 
modal  la  than  used  to  develop  an  infinite-dimen- 
sional etate  apace  nodal  for  tha  generation  of 
gyroscopic  noise. 

■semination  of  tha  literature  on  magnetic 
relaxation  (o.g.  xaf.  (91)  indicates  that  tha 
response  of  Iron  to  tranalenta  in  applied  magnetic 
field  can  be  charactarlaad  aa  tha  impulse  response 
of  a continuum  of  firat-ordor  linear  systems  with 
• uniform  volume  danaity  distribution  of  tima 
constants.  Tha  term  ‘uniform  volume  density  dis- 
tribution is  used  hara  to  mean  a spatial  distri- 
bution of  systaaw  such  that  each  volume  element 
contains  many  systems,  and  such  that  tha  systems 
in  each  volume  element  have  tima  constant*  distri- 
buted according  to  tha  same  probability  density 
function.  Each  individual  ayatam  ia  characterlsad 
by  a transfer  function  of  tha  font 

(3) 

Tha  probability  danaity  function  of  tima  constants 
(T)  la  given  by  liw  insert  in  Figure  4)  i 

I(V*n(TaA1))(VT)»  tx  < t < tJ 

0 » otherwise  I 

(4) 

Wa  shall  demonstrate  that  tha  above  density  func- 
tion la  affective  in  esplaining  tha  gyro  noiaa  ?SD 
in  addition  to  magnetic  relaxation,  which  ia  ob- 
served whan  tha  gyro  ia  operated  in  tha  presence 
of  power  aupply  transients,  incidentally,  other 
poMihla  density  forma  (for  anelastic  relaxation 
of  strain  in  crystalline  solids,  a related  pheno- 
menon) ara  dlacuaaad  by  Bowie*  and  Barry  (10) . 

The  impulse  reeponse  of  each  linear  ayatam  (Eq. 
(1))  ia  given  byi 

^(t)  - e"t/T  «> 

Tho  magmatic  relaxation  of  tha  material  is  than 
characterised  (see  Hat.  (»))  by  the  weighted  inte- 
gral of  tha  impulse  responses  of  tha  linear  sys- 
tana,  with  tha  tima  oonatant  danaity  of  Otq.  4) « 


Pd(T) 


e(t> 


■r 


(«) 


Substituting  tq.  (4) 
find  that i 

CTtT? 


Pd(T)hT(t)dT 

Bq.  (S)  Into  Bq.  (6) , wo 


m(t) 


(1/t) (a_t/T)dT 


variables , wa  obtain r 

<«  -EiT^r  IT: 


trust  the  gyro  aoins  aa  stationary.  An  explanation 


whom  wa  havn  node  tha  substitution! 
y • tA 


(7) 


(B) 


(B) 
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riMUr*  wo  obtaiai 

mW  " itt<T^T1rt,i<t/Ti> " Bi(t/tin  tlw 

Mnn  I.  (a)  la  the  eivonentlal  Integral,  defined 
byi 


»!<«»  " <V,dU 


U1J 


We  ghooM  K to  normalise  o(t)  to  ♦ (t)  > <*•«  m 
require  for  normalisation  that i 

♦(0)  - It  ♦(•)  - 0 (M) 

We  find  that i 

K - i Cl« 

Thus,  the  magnetic  dlsaccounodetioa  (raiaaation) 
ia  normalised  toi 

*(t>  - 1U<T|/t^)l,l(t/T2)',l<t/tl)l  a4) 

of  9(t),  for  t«  • 0.01,  t«  • 1.0,  on  linear- 
liaaar,  senllog,  and  log-log  aealaa  ara  found  ia 
rlgaraa  4,  5,  and  4,  respectively.  M diaowaaad 
in  naf.  (7),  *(t).  with  proper  cholaa  of  T.  and  Tj, 
oftan  fits  tha  tiaa  record  of  gyro  output  in  tha 
presence  of  power  supply  tranaiants.  Gyro  output 
ia  tha  record  of  the  torques  applied  by  tha  teg- 
aotie  gyro  torguar  in  order  to  kaap  tha  gyro  float 
Mtgle  cloaa  to  aaro.  ror  (t)  batwaan  t,  and  Tj, 
p(t)  ia  proportional  to  (-ln(t)) , a f sod  liar  raault 
la  the  a tody  of  nagnatic  raiaaation  (oaa  haf.  (11)1. 
( tocidentslly,  Tj  and  Ta  nay  ha  astinatad  by  ob- 
serving tha  gyro  output  and  using  an  analytic 
appooalaatlon  (12)  foe  bit),  for  t between  Tx  and 
t..)  to  sunnary,  tha  tins  constant  density  gluon 
la  *q-  <*>  cm  ba  uaad  to  explain  tha  deterministic 
gyro  raapanaa  to  transients.  Tha  reader  shoal I ba 
sums  that  va  do  not  haua  anplrlcal  ounfimstion 
that  tha  raiaaation  otelblted  by  gyro  output  la  tha 
prasanoa  of  power  supply  tranaiants  is  naoassarily 
nagnatic  la  origin.  Ma  can  only  ouggaot  this  as  a 
possible  source,  and  note  that  this  nadianlan  is 
effective  in  aaplaining  tha  obaaruud  power  spectral 
charaotaristlos  of  gyroooepic  output  noise,  which 
ua  shall  now  discuss. 

If  a linear  ays  tan  {eq.  ( J> ) with  tins  con- 
stant <T)  is  tod  by  an  input  function  w(t,  t)  than 
its  respond,  s(t,  t),  is  characterised  byi 


US) 


Jj  S(T,  t)  ■ t)  ♦ W(T,  t) 


Let  tha  input  function  of  tun  variablee,  w(T,  t), 
to  tha  aranaa  bo  characterised  by  oocnrianoai 
«(•)  is  tha  Dirac  dolts  function) 

*{w(t,  t)u(y,  t-ol)  4 wf(T-Y)<(o)  «•> 


w(t,  t)  is  f anally  a * tec  dtnsnslnasl  whits 
Miser,  tha  inputs  to  two  cyst  ana  with  tin  con- 
stants « and  y ara  independent  if  T F Y.  M*  HD 
My  ba  regarded  as  a stats  equation,  Mara  state 
■ft,  t)  is  a function  of  T I Hj.tjl.  Oyrosoapic 
noise  is  new  Mdslad  as  tha  weighted  integral  of 
tha  outputs  of  tbs  filters  (whom  s(T.  t)  ia  tha 
output  (at  tiaa  t)  of  a filter  with  tins 
T),  and  la  given  byi 


g(t>  - f «(t,  t)pd(t)dt 

• Jo 


(17) 


to  no  re  rigorous  form.  Kg.  (IS) , (16)  and 
(17)  are  shorthand  fori- 


g<t) 


-Jx  »d‘T>* 

./>/* 

J1,  Jo 


d|l(T,0) 


Pd(T)o"<t"*,/Td»(T.S) 


(!•) 


where  tha  first  integral,  the  initial  condition 
propagation  p ia  a wiener  integral  and  the  second 
ia  a ’ two-dlaans Iona  1 Wiener  integral*,  defined  in 
hppandla  h.  to  this  appendix  we  also  discuss  tha 
two-diasnsional  Wiener  process  0(t,s)  whose  (for- 
anl)  nixed  double  partial  derivative  is  tha  two- 
diasnsional  white  nolsa,  w(T,s) , in  eg.  (13). 
Further,  as  discussed  in  Appendix  A,  nornalisation 
of  g(t)  so  that  tha  noise  has  unit  variance  re- 
quires ■ 

(1«) 


W - lln(TJ/Tl) 


Tha  power  spectral 


ity  of  tha  nolsa  ia  then 


gluon  byi 


•gg<w  “ ’ 

11Wf(Tj  - Tj) 

1 ♦ 4«*f^TaTj. 


(20) 


Of  6gg(f)  la  plotted  lx  Figaro  7 
that  tbs  (Vf) 


A gra 


Dote 

characteristic  of  gyro  noise  ob- 
isan tally  is  iaharent  in  tha  linear- 
ised version  of  this  plot.  Tha  If2)  section  of 
Figure  1,  tha  eaperlaantally  observed  gyro  noise, 
due  to  quantisation  dominates  over  the  (1/f2)  line 
of  Figaro  7 at  high  frequencies,  seeking  .thet  por- 
tion of  tha  gyro  noise.  Further,  it  is  .fait  that 
tha  low  frequency  breakpoint  of  Figure  7 corres- 
pends  to  tines  longer  than  tha  record  lengths  cor- 
nel ly  oaployad  for  observations  of  gyro  output, 
accounting  for  its  absence  from  Figure  3 (sea  cap- 
tion of  Figure  2) . ongoing  super inents  at  The 
Charles  Stark  Draper  Laboratory  with  long  record 
lengths  of  gyro  output  indicate  that  the  pour 
spectral  density  ia  flat  at  very  lwr  frequencies 
(-1  eyole/honth)  for  sons  of  tha  gyroscopes  being 


■anna forth,  we  shall  use  tha  tern  gyroscopic 
to  refer  to  the  stochastic  process  generated 
by  o or  state  space  nodal.  Tba  gyroscopic  noise  is 
to  have  started  at  It—) , hence  to  bo 

(t-0). 

alternative  nodela,  in  terns 
of  diffusion  no  chan tana,  for  stochastic  processes 
with  tha  power  spectral  characteristics  of  gyro- 
ewaoed  ia  haf.  (11)  and  (14). 
•del  .with  a (f”l“2t)  power  spectrw,  is 
in  Mf.  (191.  Dote  however  that  because 
Miss  filter  is  a liaaar  as  tiaa  tar  only 
sated  ardsr  properties  of  tha  gyros oopic 
lafluanoa  the  tutor 


.'ll 


i i • 
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Thus  ell  mathematical  models  which  gen- 
erate atochaatlc  processes  with  tha  earns  MB  aa 
fyro  aolaa  (hanca  tha  aaaa  second-order  propartlaa 
aa  gyro  nolaa)  will  ylald  tha  awe  optlaal  (eini- 
aa  varlanca)  filter.  In  tha  neat  aactlon  wa  for- 
■ulata  tha  problem  and  indicate  how  It  can  be 
solved  In  an  infinite  diaanatonai  content. 

1.  Filtering  and  fropertiea  of  tha  filter 

Taka  X.U  real  reparable  Hilbert  apacaa. 

(0,  5,  w)  a complete  probability  apace. 


1.1 


ortho nor- 


•aparabla  Hilbert  spaea-valuod 
wariablea 


Tha  reader  la  referred  to  (20)  for  no re  de- 
tailed expoaltion  of  thla  aatarial. 

xiQ  ♦ x la  called  aa  X-valued  randoa  war labia 
(r.v.)  if  it  la  a (weakly)  measurable  nap.  Tha 
linear  apace  of  X-velued  r.v.' a la  denoted 
Met  0,  hi  X). 

Aa  X -valued  atochaatlc  procure  la  a wap 
»(.)!*♦  -»  Ma(  &.  Hi  X)~  «<•>  la  ■ aaaaurabla  pro- 

cead  If  tha  nap  (t,w)  * a(t,w)  la  awaaurabla  w.r.t. 
H,  i h (Hi  denotea  Labe  ague  naaaure  on  *♦) . 

1 a C Mat  Q,  Hi  «)  la  flrat  order  if 
n C Ll(  0.  Hi  » and  aaoond  order  if  n c L2(0,hi  *1  ■ 
for  a flrat  order  r.v.  a(w)  we  define  the  aean 
■(aim)). 


K{x(w>> 


no  r 

-»»-  h 


xtuMp 


(Bochaer  Integral) 


for  a aaoond  order  r.w.  a(w) , (h,  h)  » B(<n(w)  - 
B(ntai) ).  h>  <a(u>  - H{n(a»)  )K»)  la  a coatiauowa. 
ayawatrlc  bilinear  for*  which  haa  unique  rapreaan- 
tetlon  through  g c t(x) > Q > 0*  g*  ■ g aa  (h,K)  ♦ 
<(h,  K>.  g ia  the  covariance  of  a(u) . The  oovarl- 
of  a aaoond  order  random  variable  a(u)  la  nec- 
rily  nuclear. 

Olven  two  X- valued  aaoond  order  r.v. ‘a  ntw). 
yt«)(  (h.  K)  * B{<a(u>-x,  h>  <y(w)  - y.  K>)  haa 
unique  repreeentatlon  th , FT)  * <Xh,n>,  X c L(n).. 

H ia  called  the  covariance  of  ntw) , y(w)  and  ia 
written  oov  (ntw).  y(u)). 

ntw),  y(w)  c Mat  Q,  hi  » «•  Independent 
if  <h,  ntw)>>  A y(m)>  are  independent  for  all 
h,K  c X.  a c Lj(B,h«X)  la  oauaaian  if  <n(m),  h> 
ia  normally  diatributed  for  each  hex. 

1.2  Hlener  fioeeaa 

The  U- valued  atochaatlc  proceaa  N(t.  w)  la  a 

f lunar  oroceaa  If  11)  for  finite  colleetleme 
ti)  t RMs^]  c U,( w(tt,n),  ej)  ia  a family  of 


for  eoma  {»1>.  \ > 0 with  EA1<”« 

ami  aequance  (a.)  In  U.  WeThall  make  uae  of  the 
property  that  the  Hlener  proceaa  H(t,w)  haa  unique 
repreeentatlon 


w(t,w)  - He  £ 0 (t,w)e. 

1-1 


(limit  In  L (0.  Hi  X) 

with  tha  B.'a  Independent  real  valued  Hlener  pro- 


1.3  The  Hlener  Integral 

fuppoae  b if * - X la  locally  aaaantlally 
bounded,  maaaurabla  and  that  0(t,w)  ia  a real- 
valued  Hlener  proceaa.  Than  the  Hienar  Integral 


r 


b(t)d0(t,w) 


ia  defined  In  the  uaual  earner  aa  a limit  in 
L3(0,Hl  *1  through  a sequence  of  aimple  functions 
approximating  b(t)  In  L*(0,  Ti  X).  How  suppose 
that  b( • ) i K*  -»  l(u.  a)  satlafiea  (1)  | |a( • ) 1 1 
la  looally  eaaentially  bounded,  measurable  (li) 
t ♦ B(t)s  is  measurable  for  aaeh  a c X.  Tha  Hlener 
■tetaaral 


f 


H(t)dw(t,W) 


la  defined  In  thla  ease  aa 


■ fr 

ii-  £ 

W—  1-1  JO 


Ktle^lt^) 


real-valued  gaweeian  r.v. 'a  (11)  H(t,m)  la - 

order  for  each  t > 0 and  there  exists  soma  nuclear 
Q C t(a)  s.t. 

B{«w(t1,«d.h>  <w(tj,w)  ,h>) 

- <0t,h>  mUUj.tj) 

each  tj.tj  > 0,  h.K  c U,  (ill)  l{w(t,w))  - 0 for 

each  t > 0.  fee '((20) < P-  l*7  at  seq.)  for  proper- 
ties of*w(t,M). 

bo  ties  that  since  g ia  auelaar,  g >0, 

9*'Q 


«•) 


■ £ Vi  *V 


(limit  In  t (0,  Ti  X)) 

where  each  element  in  the  aequance  is  evaluated  as 
above,  (a.,  (t.w)  l - l,w,..  aa  msec.  3.2)  for 

•(.)  measurable  w.r.t.  tha  uniform  operator  topo- 
logy this  definition  coincides  eaaentially  with 
that  in  (120),  p.  1*0  at  seq.).  notice  that  the 
Wiener  Integral  is  defined  modulo  null-functions 
la  1^(0,  Hi  X). 

3.4  Infinite  Dimensional  formulation 
of  tha  filtering  Problem 

He  first  Show  hew  equations  (IS)  and  (1?) 
earn  be  represented  in  the  infinite-dimensional  sto- 
chastic eettlng  just  described.  Let  X - L2(t1,T,» 
R)  be  tha  spaoe  of  aquare-lntegrable  functions  with 
values  la  R,  and  let  <*,*>  denote  the  natural  sca- 
lar product  on  X.  All  random  variables  are  eon- 
sldred  with  respect  to  soma  fined  complete  proba- 
bility spaoe  (Q,  A,  f). 

Demote  by  0(t)  - w(*,t)  tha  X-valued  Hlener 
prooasa  wlth  oo variance  operator  V obtained  from 
tha  two  dimensional  Hlener  process  w(T,t). 

Consider  AiX  * Xia(T.t)  * y a(T,t) . This 

la  dearly  a bounded  linear  operator. 

bat  y(t)  • a(*,t)  be  aa  X-valued  random  var- 
glven  mo  the  solution  of  tha  integral  aqua-  _ 


* •.*»»(... 


- .•  |*il  »•  • r»* 


«Mn  To  la  •»  GwnIw  randan  wl*i« 

with  aero  naan  and  covariance  operator  1,  It  la 
aaeuned  that  y0  and  It (t)  ara  independent.  It  can 
ha  ehowu  that  (i)  y<t)  la  a Sauaeien  randan  »arl- 
ah) a.  (11)  «(y(t))  - 0 and  (111)  writlnq  A(t.a)  - 

aevlp(t).y(e>l.  <A(t,e)h.h>  - «eAtte*  *h.K*  ♦ 

r(t.a) 

«.**«•*•  W(W,,fcdo  * h.S  « l. 


»•  laanm  tha  adjoint  of  A. 

«ha  output  aquation  of  tha  yyronolea  nodal  la 

i fa 

«(t)  - I P4(T)a(T,t)dT  (**> 

Moca  fu(*)  la  a haandad  naaaorabla  function,  tha 
Moca  Mflnea  a fcaunlad  llnaar  operator  Cil  * X. 

(11)  and  (11)  oonatltuta  tha  Infinite  dinanalonal 
rapraaantatlan  of  tha  win  nodal.  Xn  practice,  tha 
fT>o  output  la  aMplad.  than,  do  lay  a ‘aaplad- 
data*  appcoo  loot ion  to  (11)  and  (11)  m ohtaln  tha 
dlaeiata-tlna  rapraaantatlon 
| ?(*♦!)  - ffM  ♦ 1(a) 


<T(n*i;  • TT<ai  ♦ bum 
»(»)  • T0 

«(n)  - jT  p^tlult.nldt  , 


n • Oilfliiti  . 

4 

■oca  ♦ c la  tha  napplap  f(t)  - a f(T) 

ahara  T > 0 la  tha  cMp  la-tine  insnaant,  and  ■(■) 
in  an  iodepeodeot  aa«nanoa  of  X-mlnad  Oanaalan 
randan  variabloo  with  aaaarlanea  oparator  9 (Mioh 
ean  ha  oalaalatad  f ran  9 and  tha  aaapllap  data) . 


* lovlata  flalta  dlMoalanal  diaorota-tlM 
llnaar  ataahaatio  oyatan  ha  fiaon  hy 

a(n»l)  - a (a)  ♦ hn(n)  Uhl 

px(n)  • h'a(n) 

■ara  e(0)  la  a damalan  Revalued  randan  aarlahla 
with  naan  0 and  eoanrlanoa  f,  a (a)  in  a "white" 

Qi#  d and  ■ ara  natrloM  of  appropclata  alaa  and  h 


a(n)  - Pj(n)  ♦ «(n)  ♦ v(n) , 
a w(n)  la  a ahlta  poaaolan  aoalar 


p,  (a)  la  tha'  ldaal  yyaa  nolaa  output  which 
■node  to  ha  aatlnatad.  in  ardar  to  aatlnata  it  wa 


haaa  to  aatlnata  a(n)  and  y(n) . thia  filtarlny 
pcahlan  am  now  ha  aolaad  weioq  otaadacd  |oflnlto- 
dlnanalonal  filtarlnp  nathada  (an,  for  aaanple, 
Pll).  ly  duality  aryunanta,  it  oan  ha  ahum  that 
thia  yaohlan  la  eqaivalMt  ta  tha  followlnp  apt  1m) 


control  pcahlan  (In  hachuard  tina) 
a(ttl)  ■ ri(t)  ♦ ou(t) > t - 

(id) 

with  ooat  functional  of  -tha  fata 

(•>(T)  ,5(T)>  ♦ £ f I p(T)«(t.t)adT 

t-0  Jlx 

♦ a(t>'Qa(t)  ♦ ou(t)2]  (Id) 

((•,•)  donotaa  tha  natural  aoalar  produoa  on 

I « R“). 

In  tha  above, 

n(t)  - (y(t).a(t))  « I « R* 
a(t)  « R 

Pit  ■ R“  ♦ I « R"  la  tha  hounded  linear  aapplny 

(t(t),  a)  » (a  T f(t),  da),  t > 0, 

•iR  *1  » I*  1*  the  bounded  llnaar  aappinq 
a ♦ (p(Ou.hu)  whom  p(-)  la  hounded  naaaurabla 
and  h I h,  I If  a poaitlwa  oparator,  0 la  a ayn- 
natrio  poaltlaa  a an  I da  finite  aatrlu,  q(x)  >e>0 
la  a hounded  naaaurahla  function  and  a > 0 la  a 
oaalar.  that  ia  of  intaraot  la  tha  aayaptotic 
hahatlar  of  thia  control  paohlon  (that  ia,  tha 
aaynptotlc  hahawlor  of  tha  filter) . 

•a  oap  that  (3d)  la  reachable  if  thore  au- 
lata  an  lntapar  y > 0 and  a oonotant  0 < a < • 


(a,  £7*0  e-r*1*)  > a| |«| |*  w * « * " R* 
!■« 


hot  m chock  whether  thia  ia  poaeibla  for  tha  I 
part  of  tha  oyatan.  No  would  require 

£ { jf  * p(T)a  7f(t)dt}J 

> a f(t)*dt  for  bum  a > 0 

and  f f t t*(T,,ti  R)  which  la  clearly  iapoaalbla. 
The  l part  af  the  oyatan  io  however  a table.  Maw 


Jf*  q(T)n(T,t)*dt  ■ «/q(*>  ■(•  ,t)  ,^F)n(*t)>. 

■once  tha  napplrp  c(-.t)  * ^P)u(*,t>  - Qu(*,t> 

nap  ha  thaupht  of  aa  an  ehoarvatloa  aquation  for 
tha  I- part  of  tha  oyatan.  *a  oap  that  the  x-part 
af  (3d)  with  the  above  ah  a creation  aquation  ia  ob- 
oaruuhla  if  * lntapar  a > 0 and  a ac actant 
d « h < • each  that 

(*.  £ r*4®  r1*)  > h| Ul I1  *m. 

lot 

f Inca  q(t)  > c » • tha  l-part  of  the  oyatan 
io  ahaaraahla.  now  aoooa  that  a(k*l)  - da(k)  ♦ 


io  otahlllaahla  and  £d,k®dk  > 0 for  aoho  o. 


t*4  l 


AmIIhm  I'l.  .«■  (Hit  |..er  nuudKi  of  y.uir  | ■ <'ci. 

I XI 


; *» . 


Xt  than  follow*  from  th*  reeulta  of  Hager  and  Horo- 
wita  (22 1,  on  th*  aayuptotlc  bahavlor  of  Dlacrata 
Riccati  Operator  Equation* , that  tha  reaulting  fil- 
ter la  aayeptotically  atabla.  ror  tha  appropriate 
oonoapta  of  filter  atablllty  aaa  tha  forthcoming 
by  Vinter  (2)1. 


Appendix  A 

Tha  Two-Dlmenaional  Wiener  Procaaa 

In  thla  appendix,  we  diacuaa  the  two-diann- 
elonal  Wiener  proeeae.  with  cover laneai 

E(0(T,t)0(e,o))  -W  • mind. a)  • mln(t.o) 

(A-l) 

Mote  that,  formally,  the  mixed  double  partial  de- 
rivative of  thla  procaaa  will  have  the  covariance 
•f  * two-dlmaaaional  white  noiae.beeauae,  formally! 

a3  a2 

BUTOt#(,'tm?3S*  #u»°m 

»4 

• E(0(T,t)0(e,o) 1 (A-2) 

and  from  tq.  U-l)  we  have  thati 

.4 

>t>t?oS  *C0(T.t)0(a.o>l  - w«4(t— a)*d(t-o) 

(hr)) 

We  ehall  firat  ahow  th*  exletenc*  of  a (Gauealan) 
random  procaaa  with  th*  covariance  in  Eq.  (A-l) . 

We  than  ahall  give  meaning  to  a two-dlmenalonal 
Wiener  Integral  (where  f C Ll((0,*)  * (0.  •)))• 


jfr 


f (T.t)dS(T.t) 


(A-4) 


rinally,  we  ehall  ahow  that  th*  model  w*  have  given 
for  gyroecoplc  no lee  in  Sq.  (Id) ■ 


•<t> 


- jX2p4(x)«~t/1av<r,oi 


*n‘v 


(T)*"(t‘">Ad6(T.n) 


IA-S) 


in  fact  yield*  the  dealred  power  epectral  denalty 
(■q.  20)i 

V"  - ‘imftTqy  I5F  * 


t 


,l|2vfd2  - V 


«A-4)  * 


Further  detail*  concerning  th*  nultlp* matter  Wie- 
ner prooeee  and  Wiener  integral  are  diecueeed  by 

( 1. 

WO  firot  ahow  th*  emietenoa  of  th*  two-dlnan- 
lanor  prooaaa.  The  argument  oloaely  par- 
allel a that  of  J.N.C.  darn  (17) . Cheoae  two  aet 
•f  oanpleta  erthonocnal  fueotiona  in  lJ<(0,  -)). 
re  <•»•>  la  aoalar  pro  duct  notation) 

(♦4<t»>  . . <dt.  Y • j£"  ♦1(t)fJ(t)dt  - 


(o » vj) 

■ ( } i 1-1 . 2.  ...  IA* 

U»  1-JJ 

{♦4(t>)  . . <d4.  ij*  - 7)0j(7)dt  - 

{°l  1**J> 

* f i*&*  2*  •••  ca- 

ll i-j» 

(Wot*  that  (♦,(!)}  and  (*.(T))  may  be  the  earn* 
eat  of  function*.)  Next,  define  a aequanc  of 
doubly- Indexed  Gauealan  random  variableai 


t*lj>  • • BUijV  " VV 


(1-1,  2,  ...\ 

J— 1.  2,  .«./ 


(A -9) 


where  4a_  in  tha  Kroneckar  delta  function,  defined 

»»y 

- r*  " j <A-io) 


i°, 

111  i— J 


We  now  define  a aequanc*  of  random  proceeaee 
(0*(T,t)}.  Each  random  procaaa  la  a function  of 
two  variable* , (t)  and  (t) . The  definition  la 
given  byi 

0^(T,t)  - J ♦1(Y)f<(a)dndY 

iJ°J°  1 (A-ll) 

«x  (T)  and  (t).  W*  claim  that  {*"(7,0))  la  a 
quadratic  naan  Cauchy  convergent  aequanc*  of  ran- 
dom variable*.  Obeerve  thati  (aay  M > W) 


■ ((^(t.t)  - aNt.t))2! 


(a)dodY  * 


u.a-.jccv- 


12) 


*y  Eq.  (A-*) , we  obtain i 

W((0"(T,t)-«*(T.t))a)  - W. 

[ £ ( * ♦1(Y)«Y)al*[  £ < 

[limn  0 1 J 11-4K1  ■ 

Define  th*  following  function  i 

(A-12) 

lit  *«Y<t| 

xtm  - < I 

* (Oi  Y > * * 

(A-14) 

■q.  (A-l)  in  dot  product  nota- 
*C(l^(t,t)-«"(T.t))a)  - 

.1 


(A-iJ) 


...U 


•!WH  I i " (Kll 


: i 

t 


J 


By  tlM  ettlmwcMl ity  of  the  sequences  (0.  (Y>>  and 
(O.(X))  we  have  that  aach  factor  on  tha  right  of 
Iq.  (A-1S)  approach aa  taro  aa  M ♦ -.  Thus 
{f*(X,t))  ia  a quadratic  Man  Cauchy  convergent 
a a ni>a  nr  a of  randan  variables.  Call  tha  Unit 
8{X,t).  Wa  nou  daaonatrate  that  8(X.t>  haa  oovari 
anco  aa  in  Iq.  (i-ll • . 

It$(x,t)0(s,o)1  - «-|£<VIt><W| 


f(X,t)dB(T,t)-  £al  I i 

vi  yoyo 

/ |f  (T,t)-f(T,t)|2dTdt 
i Ia  ” 


f (T,t)dfl(T,t) 


f (T,t)d0(T,t)afin  in  q.n, 


B(6(x,t)6(s,o))  ■ W*nln(t,a)*ain(x,s>  U-17) 

No  have  thua  eetabl  iahad  tha  existence  of  tha  Me- 
diae nsional  Nlanar  procaaa. 

Na  now  wish  to  define  a two— rt loan  a tonal  Wiener 
Integral.  We  firat  develop  an  analog  of  the 
'orthogonal  lncranent*  property  of  Brownian  Motion. 
Mr  a Brownian  notion.  v(t),  it  ia  well  known  (1?) 
that  i 

U(-)  denotea  length) 

Bt(li(t)-V(e))  (V(q)-M(r) ) )■*•!( (a, t)  Olr.q)) 

U-18) 

Define  the  following  function,  over  a box  (lx2,x2) 

■ (tL,t2))  in  the  firat  quadrant  (no  generality  ie 
loet  by  defining  our  procoaa  for  x > 0,  t > 0)  of 
the  K2  planai 

r(Tl.T2.tl.t2)-e(Ta.tJ)-8(t2,tl)-8(Tx.ta) 

♦ Ktj.t^  W- it) 

Bote,  incidentally,  that  if  a deteminiatlc  func- 
tion q(T,t)  c CJ  had  ite  nixed  do«d»la  partial  de- 
ravativa  integrated  over  thle  donaln,  we  would 


• (X,t)d0(X,t)  <h-23! 

The  claae  of  fvaictiona  f(T,t)  for  which  thla 
ia  poeeible  ia  1.2(10,-) *10,-)) . In  addition,  aa 
in  the  cno-dlMne lonal  case  in  Wong  (IB) , wa  find 
thati 


g(T,t)d£(T,t) 


f(x,t)d6(x,t) 


f(x,t) (g(x.t)dxdt 


We  shall  Make  use  of  Cq.  (A-24)  in  showing 
that  our  nodal  for  gyroscopic  noise  (Eq.  (A-S) 
yields  the  desired  power  spectral  density  (Bq. 
(A-d) ) . fron  Bq.  (A-5)  wa  have  that  (for  o > 0) i 


*lg(t)g(t-o)J 


(t-oj/l 


q(x1.«2>  ♦ q(xx.tx)  (A-20) 

in  be  taken  as  activation  for  Iq.  (A- 1*1 . It 
Lly  aeon,  using  only  Iq*  (A-17)  that!  (A(>) 


|(t-W)-Bl 


• W*  At ( txx ,Xj) * (tx ,tj) ) O (tXj.X4)*ttj,t4))J 

(A- 21) 

ton  (A- 21)  to  the  analog  fee  Iq*  (A-ll).  doing 
ereeerty,  we  paweeed,  as  Wong  (1?)  dees  for  a 


g(y)dy  (Y,0) 


f (x)d(i(X  ,0) 


f (x)g(x)  (h(x)dx) 
((or 


Mm,  m haw  that  0(T,s)  (T  > 0,  a >0)  la  inde- 
pendent of  u(t,0)(t  > 0).  Thus  w*  obtain  fro* 

Bq.  (A-3S)  that i 

B(g(t)g(t-<*)1 


h(T)dT 


♦ pl<v'2<t’a,/Tfo 

Integrate  In  Bq.  (*-27)  to  obtain 
B(«(t)9(t-a)) 

" * 4(T)."(t^,/T  h(T,dx 

\ 

a .^JH*  ,«)p*(T)U^-*<t^>AJ4l 


s2,/T  MdadT 


(*-37) 


(*-38) 


Lot  t * • in  Bq.  (*-28)  .to  obtain  a stationary  ran- 
do*  procasa  charactarixad  byi 


B («)  - E(g(t)g(t-al) 
99 


- f2  ,2t,.-l“lA  pj|(x)dr 


(*-38) 


* Fourier  transom  of  Bq.  (*-29)  yialdai 


S <f)  - 


f (-tttK 

Jo  \WiTT  / 


(T)PdT  (*-30) 


itutiag  Bq.  (4)  Into  Bq.  (A- 30)  and  using  a 
lxatlon  (M  • 2h(t]/i1l)  to  give  8 (f)  unit 

ca,  wa  obtaini  w 
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